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Abstract 



The transverse momentum dependence of J/ tp suppression in rela- 
tivistic heavy ion colhsions is calculated from initial state gluon rescat- 
" . tering not only with nucleons but also with prompt gluons produced in 

. nucleon-nucleon collisions in the early phase of the reaction. The sec- 

. ond mechanism is new and turns out to be important in a description 

. of the data for central collisions. 

o 

1 Introduction 

o 

C ■ Recently, the NA50 collaboration has presented the final analysis of trans- 

verse momentum distributions for charmonia produced in Pb-Pb collisions 
X ; at 158 A GeV. For the J/"^, the values of {pt), {pi) and (M^ - M) exhibit 

^ ' a similar trend as a function of the centrality of the collisions: they first in- 

crease and then flatten when the centrality increases, and finally turn steeply 
upward for the most central collisions. 

It is commonly accepted^] that the new data on J/\I' suppression in 
Pb-Pb collisions cannot be explained only by inelastic collisions of the J /'^ 
with the nucleons of projectile and target. The anomalous suppression needs 
a new absorption mechanism, and the observed drop at large Et (~ 100 
GeV) is either considered as the expected onset of J/\& melting]^, Q in the 
picture of the quark-gluon plasma, or interpreted H, ^ as an effect of 
the transverse energy fluctuations around its mean value determined by the 
collision geometry. The data on J/\I' transverse momentum distributions 
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provide us additional interesting information and possibly a deeper insight 
into the underlying physics. 

In a recent publication Q we have drawn the attention to a close cor- 
relation of the values of {pf){Et) and the degree of anomalous suppression 
AS{Et) = SG{Et) — Sob{Et), where SciEt) (G= "Glauber" ) is the suppression 
calculated from the inelastic scattering of the J/ip on nucleons alone, while 
"06" refers to the observed values. In the same paper we have proposed one 
explanation for this correlation: The escape of high pt J/ ip from final state 
absorption in a QGP or an environment of comovers. Unfortunately, the the- 
ory required a new parameter. In the present paper we take an alternative 
point of view: We investigate whether the data on {pf){Et) can be described 
by an anomalous mechanism in the initial state . We have followed this line, 
since the bulk part of the observed (pf) seems to be generated in the initial 
state. 

It is generally agreed that the additional transverse momentum broad- 
ening Sp"^ = {pI)ab — {vDnn observed in nucleus-nucleus (AB) collisions 
over the one in nucleon-nucleon (NN) collisions arises in the initial state, i.e. 
before the production of the cc pair[|T^. The two partons (mostly gluons) 
traverse nuclear matter thereby colliding with nucleons (thus acquiring ad- 
ditional p1) before they fuse to the cc. Indeed, the data for {p1)ab show a 
linear relation with the path length of the gluons in the initial state as is well 
studied in pA and S — U collisions[||. In a nucleus-nucleus collision, espe- 
cially with heavy nuclei, an additional phenomenon arises, the phenomenon of 
wounded nucleons. In the approach of prompt gluons proposed in refs. [^, |3 
for the description of anomalous J/-?/' suppression "wounded" nucleons are 
nucleons which contain additional prompt bremsstrahl gluons whose number 
is proportional to the number of inelastic A^A^ collisions. In the approach of 
refs. [1^ the additional gluons lead to a larger inelastic cross section for a 
J/ip on a. wounded nucleon compared to a ground-state nucleon. 

The wounded nucleon approach should also have an effect on the initial 
state. The gluons from projectile and target nuclei, which eventually fuse 
to form the cc, encounter nucleons and wounded nucleons before they fuse. 
The interaction of these gluons with nucleons is identified as the main source 
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of the additional transverse momentum broadening 6pf. If the encountered 
nucleon is wounded, 6p^ should be even larger. It is this effect, which we 
calculate in this paper. 

In our calculation we rely heavily on the notation and on the parameters 
chosen in refs.[0, yet our presentation is essentially self-contained. We 
first present the formulae for J/ip suppression in Pb-Pb collisions calculated 
from inelastic J /ip N interactions and from collisions with prompt gluons in 
the final state. Then we turn to the initial state and derive the formulae for 
the generation of the transverse momentum (Sect. 2). Then, in Sect. 3, we 
present the numerical results and compare with experiment, followed by the 
conclusions (Sect. 4). 



2 Generation of additional transverse momen- 
tum via interaction with prompt gluons in 
the initial state 

Fig. 1 is a schematic representation of J/il) production in nucleus-nucleus 
collisions with particular emphasis to the prompt gluons in the initial and 
final state. The collision of a row of projectile nucleons on a row of target nu- 
cleons is depicted in a two-dimensional (time-longitudinal coordinate) plane 
and in the A^A^ c.m.s. A gluon from a projectile nucleon and a gluon from 
a target nucleon fuse to produce a premeson at point O which has three di- 
mensional space coordinates (s, za) in the rest frame of the projectile A and 
(6 — s, Zb) in the rest frame of the target B, where h is the impact param- 
eter for the AB collision. The premeson propagates through the projectile 
and the target, and develops in time towards the asymptotically observed 
J/\E'. During the propagation the premeson interacts with nucleons from the 
projectile and target. Since these nucleons have experienced collisions with 
other nucleons before their encounter with the cc premeson, and since they 
have radiated bremsstrahl "prompt" gluons (dashed lines), the cc premeson 
will also interact with these radiated gluons. Taking into account the two 
kinds of final state collisions of the premeson the J /'^ suppression function 
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can be written as 



d bd sdzAdzBK{b, s, za, Ei 



t) , 



where Et is the transverse energy which reflects the centrahty of the collision, 
and the expression for the suppression kernel K is taken as||lCI||: 



K{b, s, ZA, ZB, Et) = pa{s, ZA)pB{b~ s, ZB)e-^'^^'^^P{Et\b) 



(2) 




Figure 1: A two-dimensional (time-longitudinal coordinate) plot for J/\E^ pro- 
duction in a collision of nuclei A and B in the c.m.s. of the colliding nucleons. 
Nucleons are denoted by thin solid lines, the two gluons which fuse to form 
a premeson at point O are denoted by thick solid lines, while the bremsstrahl 
gluons are drawn by dashed lines. 

Here PA,B{'f^) is the density of the nuclei A and B, respectively. The final 
state J/\E' absorption via the interaction with nucleons and prompt gluons 
in the geometric row labeled by the variables b, s is contained in the two 
exponentials given by 

lA{b,S,ZA,ZB) = 
+ 

lB{b,S,ZA,ZB) = 
+ 



dzCTab^ it^)pAis,z) 



^il^ / dz dz'aJ^{t^)e{tg){ng){tg)pAis,z)pB{s-b, 



a, 



in 

-oo 
NN 



dza^fl^(tq,)pB{s- b, zb) 



ZB 



dz 



dz'oll (txt ) (tg) {ng) {tg)pB{s- 6, z) p a{s , z 



where is the inelastic NN cross section, (n^) is the mean number of 
prompt gluons produced in a NN colhsion, and af^ and cr*f,^ are the J/\l/ 
absorption cross sections for its interaction with a nucleon and with a prompt 
gluon, respectively. Effects of formation time for the development of a cc 
pair to the J/ip and for the radiation of a gluon are taken care of by the 
time dependence of cr*j,^ (t^) and {ng)(tg), respectively. Here is the time 
difference between the "^N interaction point and the formation point O of the 
premeson, and tg is the time difference between the "^g interaction point and 
the NN interaction point which produces the prompt gluon. The 6 functions 
in Ia and Ib ensure that only the prompt gluons, which are created before 
the interaction with the effective J/\E', are taken into account. 

Using the two-channel model for the evolution of a cc pair, the time 



dependent effective cross section a^^J^ can be expressed as[^ 



C'lir) = + «^ - ) cos(AMr) , (4) 

where af^^ is the initial (r — > 0) and cr*^ the asymptotic (r oo) absorption 
cross section, and AM is the mass difference between J/'^ and We use 
values cr*^ = 6.7 mb and af™ = 3 m&p]|]. 

In the constituent quark model, the inelastic "^g cross section differs from 
the inelastic "^N cross section by a factor 3/4 
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Within pQCD the time dependence of the mean number (ug) of the radi- 
ated gluons can be calculated with an adjustable parameter Umin, the soft 
limit of the transverse momentum carried by the prompt gluons. In order to 
reduce the numerical efforts, it is parameterized as a step function 



{n,)it,) = { (6) 




In our numerical calculations we have used fixed values|]TO[ to = 0.6 fm/c 



and n° = 0.75 (corresponding to Umin = 0.75 GeV) for Pb-Pb collisions at 



g 

158 A GeV. 
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To include the effect of transverse energy fluctuations which have been 
shown significant for the explanation of the sharp J/\E' suppression in the 
domain of very large Et values, we relate the mean number (ug) to the 
transverse energy in the same way as the energy density]^] and the comover 
density]^, 

K) - , (7) 

where the mean value {Et){h) is proportional to the number of participant 
nucleons at fixed 6, {Et){h) = qNp{b), and Np{b) is determined by the collision 
geometry. The function P{Et\b) in the suppression kernel (H) describes the 
distribution of transverse energy in events at a given impact parameter b. It 
is chosen as a Gaussian distribution [§] 

PiEAb) = =e ^i-'-^vW (8) 

/27rg2aiVp(6) 



with the parameters q = 0.27AGeV and a = 1.270. 

We now turn to calculate the J/"^ transverse momentum within the mech- 
anism of gluon rescattering in the initial state. Before the two gluons, one 
from the projectile and one from the target, fuse to form the cc premeson at 
point O, either of them may scatter off one or several nucleons and one or 
several prompt gluons, as shown in Fig.l. Via the exchange of transverse mo- 
mentum at each interaction vertex, the mean squared transverse momentum 
of the J/\E' is a function of the history of the two gluons. 



{pi){b,S,ZA,ZB) = {pt)NN + ^^Y^^9Nib,S,ZA,ZB)+ , f iUb, S, Za, Zb) ■ 

^gN ^NN^gg 

(9) 

Here the first term is the contribution from the cc production vertex in an 
isolated NN collision, while the second and third terms arise, respectively, 
from multiple gN and gg collisions before fusion. Since the number of the 
gN collisions is proportional to the number of nucleons which the two gluons 
encounter, the second term scales with the sum of the lengths of the two 
gluon paths, 

igNib,s,ZA,ZB) = dzpAis,z)/po+ dzpB^s - b, z) / pQ . (10) 

J —CO J zb 
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The constant before the length in eq. (|^) depends on the mean transverse 
momentum {pDgN acquired in a gN colhsion and on the mean free path 
XgN of a gluon in nuclear matter. The third term in eq.(y) is proportional 
to the number of A^A^ collisions in which the prompt gluons are generated. 
Therefore, is a sum of two 2-dimensional integrations (similar to the 
second terms in I a and Ib in the final state interaction), 

lgg{h,s,ZA,ZB) = dz dz'{ng){tg)pA{s,z)pB{s-h,z')/pQ 

J zb J —CO 
Zyl rco 

dz dz'{ng){tg)pB{s-b,z)pAi7z')/pl- (H) 

-co Jzg 



The other factors in the last term of eq.(|^) are: The mean free path A} 
Po'^NN inelastic A^A^ scatterings (the source of the prompt gluons) and 
A~g = Po<Jgg the mean free path for the interaction of a gluon (which fuses to 
form the cc) with a prompt gluon. We use aj^'^ = 32 mb and po = 0.17//m^, 
while {Pt)gg/Xgg is Considered as a free parameter. The effect of transverse 
energy fluctuations can also be included in the initial state interaction by 
making the transformation (|^) for the mean number (ug) of the prompt 
gluons, like in the final state interaction. 

We will call i'^g the squared gluon length related to gg interactions. How- 
ever, it is necessary to note that while ig^ is a geometrical length, either 
term in i'^^ is not a simply product of the two geometrical lengths, since pro- 
duction time effects of the prompt gluons are included via the mean number 

(%)• 

The expression (|^) for the mean transverse momentum holds for a J/\E' 
produced from the rows with coordinates b, s, za and zb- In order to compare 
with the experimental data, we have to average. The average of a quantity 
Q(b, s, Za, Zb) at fixed transverse energy Ef is calculated from the suppression 
kernel K (eq.(D) 

{Q){Et) = I d^bd^sdzAdzBQib, s, ZA, ZB)K{b, s, ZA, Zb, Et)/S{Et) . (12) 



Inserting expression eq.(0) into eq.(O) one obtains 



{lgN}{Et) + 

^gN ^NN-^gg 



{Plm) = {p',)NN+'-^{igNm) + i^{fgg){E, 
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{p't)NN+^-^{f/Jm), (13) 



where we have introduced the effective gluon length {i'g{J) in nuclear matter 
defined by 



{VtlgN/XgN 

3 Results and discussions 

To simplify the numerical calculations we use a constant nuclear density 
p(r) = pqQ{Ra — r). Since the calculation of the average transverse mo- 
mentum needs information of the suppression function S{Et), we first cal- 
culate the ratio of the J/\E' to Drell-Yan cross sections a^,{Et) / aDY^Et) = 
S^{Et)/ SdYiEt). Here 5"^, is the J/\E' suppression function eq.(|l|), and the 
Drell-Yan cross section Sdy is obtained by putting = = in the kernel 
K of S. The results are shown in Fig. 2 as a function of Ef together with the 
new data of the NA50 collaboration. The calculation is normalized a.t Ef = 
(NN collision) to crif,/aDY = 50. In the absence of the Gaussian distribution 
in Et (eq.dl)) and its fluctuations related to {ug) (eq.(0)), we have the same 
results as in||lO|. It is obvious that the final state interaction with only nucle- 
ons (dashed line) cannot reproduce the data. The results including final state 
interactions with prompt gluons (thin solid line) fit the i?f-dependence well, 
except for the fast drop at very large Et values. According to refs. ||, |12 



the drop at high Ef arises from transverse energy fluctuations. Including the 
Et fluctuations via eq.(0) the agreement with data is improved, but less than 
the data seem to require. The quality of the fit is about as good as in ref. 
and the residual discrepancy in the agreement with data may have its origin 
in the Et assignment (c.f. ref.pj). 

The average gluon lengths (eq. (p^) are plotted in Fig. 3 as a function 
of Ef. The length (igN) in nuclear matter (dashed line) and the length 
y {£gg) related to prompt gluon contribution but without consideration of 
transverse energy fluctuations (thin solid line) have the same behavior and 
almost the same amplitude. They increase in the domain of small transverse 
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Figure 2: The relative J/\l/ production cross section (T^,/aDY in Pb-Pb colli- 
sions at 158 A GeVl\r^ as a function of transverse energy. The dashed line 
corresponds to the normal nuclear suppression, the thin solid line is obtained 
with final state rescattering with nucleons and prompt gluons, and the thick 
solid line is calculated including the fluctuations in transverse energy. 

energies [Ef < 50 GeV), and then saturate with Ef for values Ef > 50 
GeV, where one observes large anomalous suppression and the prompt gluons 
become important. The saturation of the gluon lengths and of the transverse 
momentum (shown below) in this domain is due to the saturation of the 
amount of matter traversed by the initial gluons in the Pb-Pb collisions. 
As for the suppression in Fig. 2, the contribution of the transverse energy 
fluctuations is significant only for very high transverse energies, Et > 100 
GeV (solid line in Fig. 3). 

The effects of initial gN and gg scattering to the average transverse mo- 
mentum are controlled by the two parameters, {Pt)gN/^gN and {vDgg/^gg^ 
respectively. The former may be extracted from experiments with light 
projectiles, and the latter is determined from the fit of eq.(|) to the data 
in the Et region where prompt gluons become important. With the iso- 
lated NN contribution {p1)nn = l-H {GeV/cY and the gN contribution 
{Pt)gN/^gN = 0.081 (GeV/c)'^/ fm^, the transverse momentum calculated 
with only gN scattering (dashed line in Fig.4), namely considering only the 
first two terms in eq.(^, agrees well with the data in the domain of Et < 50 
GeV, while it is clearly below the data outside this domain. Including the 
initial scattering with prompt gluons but without the factor eq.(|^ (thin solid 
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line) and choosing the parameter {vDgg/^gg = 1/35 {Pt)gN/^gN, satisfactory 
agreement is obtained for all the data except for the last point, where the 
transverse energy fluctuations are expected to play a crucial rule. We include 
the fluctuations in Et via the factor eq.(|^) which becomes relevant only for 
i'^g in eq.(^ and arrive at the solid line in Fig. 4. There is some improve- 
ment, but the significance is difficult to judge, since there is only one data 
point with rather large error bar. 

12 
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Figure 3: The gluon lengths {igN) (dashed line) and \J {^1g) without (thin 
solid line) and with (thick solid line) fluctuations in transverse energy as a 
function of the transverse energy. 
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Figure 4: The average squared transverse momentum (pf) as a function of 
transverse energy. The data are fromj^. The dashed line corresponds to the 
normal gluon rescattering on nucleons, the thin solid line is obtained with 
gluon rescattering with nucleons and prompt gluons, and the thick solid line 
is calculated with also the fluctuations in transverse energy. 
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In order to see the results of our calculations in a different form, we display 
the experimental and calculated values for {pl){Et) as a function of {igN){Et), 
Fig. (5a), and as a function of {t/J)iEt) (eq.(|n|)). Fig. (5b). While the de- 
pendence (pf) versus Et is rather complicated with strong saturation effects 
for large values of Et, the dependence (p^) versus (igN) or (tg^J) is nearly a 
straight line. The dashed and solid lines in Figs. 5a,b correspond to the theo- 
retical curves in Fig. 4. According to the definition eq.(0), the slopes of the 
two straight lines are the same and equal to {Pt)/^gN = 0. 081 {G eV / c)'^ / fm. 
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Figure 5: The observed {pf){Et) plotted versus the calculated {igw) and {(^g^) ■ 
The data are from^. 

Until now, we have studied the behaviour of (pf) as a function of Et. In 
the experiment under consideration [0] the authors have also measured the pt 
distributions within a given Et interval. They have displayed their data as a 
ratio 



where S{pt\Ef') is the differential distribution inp^, normalized to / d'^pt S{pt\E 
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S{E^^). We have calculated S{pt\E^/^) by assuming a Gaussian pt depen- 
dence. Then the suppression function S'(pj|£'f'') can be written as 



S'(pt|£'f*^) = / d'^bd'^sdzAdzBK{b, s, za, zb, 



1 



{Pt)ib,S,ZA,ZB] 

(16) 

For the comparison with the experimental data, we divide the Ef region into 
five bins with mean values Et = 23, 48, 70, 90, 110. The theoretical results 
together with the data are shown in Fig. 6. 
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Figure 6: The ratio Ri/i of the J/\E' cross section with large value Ef in the 
Et bin i to the cross section with small e[^^ in the first Et bin as a function 
of Mt — M. We have plotted the logarithm of this quantity, since the global 
data seem to be exponentials exp(— (M^ — M)/T). There are five Et bins 
corresponding to mean Et = 23, 48, 70, 90, 110 GeV . The data arefrom^J. 
The theoretical lines are calculated with initial and final state rescattering 
with nucleons and prompt gluons and with transverse energy fluctuations. 
For -R5/1, the solid line is our calculation, while the dotted line represents the 
calculated curve normalized to the data for small Mt — M. 

For the figures we have not followed ref. H], who plot Ri/i versus pt, but 
we rather plot lni?j/i versus AM = + Pt ~ M. This representation 

leads to a straight line for the data integrated over Et (cf. |[l[]). According 
to Figs. 6, this representation also leads to nearly straight lines. The solid 
lines represent our calculations. For all ratios, except -R5/1, the data are 
well accounted for. We see a problem in the quantity -R5/1. Our prediction 
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is consistently below the data, although our integrated value S{El ') with 
Ej.^^ = 110 GeV lies above the data (Fig. 2). Thus we would have expected 
our calculation to be consistently above. In order to see whether calculated 
and experimental slopes agree we have normalized the calculated curve at 
small pt to the data. The renormalized curve is shown by the dots and de- 
scribes the data reasonably well, though there may be systematic deviations 
at large pt- But since the error bars are too large, speculations about the 
discrepancies may be dangerous. 



4 Conclusions 

In this paper, we have studied the transverse momentum dependence of J/ ip 
production in Pb-Pb collisions at 158 AGeV. We have investigated whether 
the mechanism of prompt gluons can account for both effects, anomalous 
suppression of J/ip and anomalous values of {pf). Here, "anomalous" is 
judged with respect to the "conventional" values of suppression and {pf) 
derived from the data of pA and S — U collisions. The conventional physics 
is shown as dashed lines in Figs. 2 and 4, while the anomaly is the difference 
between the data and these lines. While the anomaly is large (an effect of a 
factor 2) in the suppression data. Fig. 2, the anomaly is small for the values 
of {p^). The difference between the data and the dashed line is of order 10%, 
if we use {pI)ab — {vDnn as the baseline. Within the model of prompt gluons 
we are able to describe the data for the J / ip suppression (except possibly at 
very high values of Et) by introducing one adjustable parameter 

. J/4'g 

^77^ = 0-56 (17) 

'^abs 

where is the number of produced hard prompt gluons per A^A^ collision 
(eq.(||)) and the cross sections (^ai^^ and o'^^^ describe the inelastic cross 
sections leading to the destruction of the J/ip in collisions with a gluon and 
a nucleon, respectively. Within the same model we are also able to describe 
the data for (pf), Fig. 4, after a new parameter 



1/35 (1^ 
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is adjusted. The magnitude of the values eqs.(0) and (p!8| ) for the two 
parameters reflects the observation that the anomaly in the suppression is 
large (50%) while it is a small effect for the values of {pi). 

We conclude: We have presented two possible explanations for the anoma- 
lous behaviour of the dependence of {p1){Et) in Pb-Pb collisions. In ref.|]^ 
the effect is described as one of the final state while in this paper the gluon 
rescattering on bremsstrahl gluons in the initial state has been investigated 
as a possible source. In both approaches, satisfactory agreement with the 
data is achieved after adjusting one parameter in each calculation. 
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